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ABSTRACT 



We consider the tidal excitation of modes in a binary system of arbitrary eccentricity. 
For a circular orbit, the modes generally undergo forced oscillation with a period equal to 
the orbital period (T). For an eccentric orbit, the amplitude of each tidally excited mode 
can be written approximately as the sum of an oscillatory term that varies sinusoidally 
with the mode frequency and a 'static' term that follows the time dependence of the tidal 
forcing function. The oscillatory term falls off exponentially with increasing b a (defined as 
the ratio of the periastron passage time to the mode period), whereas the 'static' term is 
independent of b a . For small b a modes ib a 1), the two terms are comparable, and the 
magnitude of the mode amplitude is nearly constant over the orbit. For large b a modes 
ib a > a few), the oscillatory term is very small compared to the 'static' term, in which case 
the mode amplitude, like the tidal force, varies as the distance cubed. For main sequence 
stars, p, /, and low order (/-modes generally have large b a and hence small amplitudes 
of oscillation. High overtone (/-modes, however, have small overlap with the tidal forcing 
function. Thus, we expect an intermediate overtone (/-mode with b a ~ 1 to have the 
largest oscillation amplitude. In addition, we find that the mode amplitude is independent 
of the dissipation rate except when the mode frequency is very close to orbital resonance or 
the damping time is less than T; both conditions are unlikely. Moreover, orbital evolution 
causes a resonant mode to move off resonance with time. This severely limits the amplitude 
of modes near resonance. Rotation of the star shifts the mode frequenices but has otherwise 
little effect on the mode amplitude (provided that the rotation rate is small). Hence, tidally 
excited modes have amplitudes and phases that are periodic with period T, making them 
readily distinguishable from oscillations excited by other mechanisms. 

We apply our work to the SMC radio pulsar PSR J0045-7319, which is believed to 
be in highly eccentric orbit with a 10 Mq B-star. We find that the g-j mode (with period 
of 1.07 days) of the B-star has the largest oscillation amplitude, with a flux variation of 
2.3 m-mag and a surface velocity of 70 m s _1 . The flux variation at periastron, summed 
over all modes, is about 10 m-mag; in addition, we propose that the shape of the light 
curve can be utilized to determine the orbital inclination angle. The apsidal motion of this 
system, calculated without the usual static approximation, is larger than that predicted 
by the classical apsidal formula by about 1.0%. For the PSR B1259-63 system, the tidal 
amplitude of the Be-star companion is smaller by a factor of 70 due to its larger periastron 
distance. To understand the dependence of tidal excitation on stellar structure, detailed 
numerical calculations of modes of a general polytropic star are also presented. 
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§1. Introduction 

The study of tidal interactions in binary systems is, among other things, useful as 
a test of stellar structure (the apsidal motion test, cf. Schwarzschild 1958; Claret and 
Gimenez 1993), as a mechanism for the formation of close binaries in globular clusters 
(Fabian et al. 1975; Press and Teukolsky 1977; Lee and Ostriker 1986; Kochanek 1992), 
and for explaining the synchronization and circularization of orbits (Lecar et al. 1976; 
Zahn 1977, 1989; Hut 1981). The primary focus of this paper, however, is the study of the 
tidally induced oscillations of a star and their use as a possible probe of stellar or orbital 
parameters. 

Our work is in part motivated by the recent discovery of two radio pulsars in binary 
systems with main sequence stars (PSR J0045-7319, see Kaspi et al. 1994; and PSR B1259- 
63, see Johnston et al. 1994). In the past, radio pulsar timing observations have been used 
to determine orbital parameters with unprecedented precision (Taylor and Weisberg 1989; 
Kaspi et al. 1994) Since both stars in the binary systems were compact, tidal interactions 
could be neglected. Indeed, the absence of tidal effects in binary pulsar systems has 
allowed various general relativistic effects to be measured. The discovery of binary radio 
pulsars with normal star companions, however, has opened up the possibility for measuring 
tidal interactions with great precision to probe the interior of normal stars. PSR J0045- 
7319, because of its high eccentricity and small periastron distance, is currently the best 
candidate for finding tidally excited oscillations. In fact, PSR J0045-7319 has one of the 
smallest values for the periastron distance (in units of stellar radius) of all the known 
detached binary systems. 

Tidal deformation of a star can be modeled as a linear superposition of its non-radial 
pressure and gravity normal mode oscillations (Press and Teukolsky 1977). We investigate 
how individual mode amplitudes vary with modal and orbital parameters for a Keplerian 
orbit with arbitrary eccentricity. The dependence of the mode amplitude on the rotation of 
the star will also be discussed. We ignore, however, the new family of modes that appears 
in a rotating star (see Papaloizou and Pringle 1978, for a discussion of the rotation modes). 

We calculate the surface oscillation amplitude of the star as well as the perturbation 
to the pulsar orbit as a possible means for detecting these oscillations. We also calculate 
the detailed shape of the light curve near periastron, which can be used to determine the 
orbital inclination angle for eccentric binaries such as PSR J0045-7319. 

The plan of this paper is as follows. In the next section we establish the main equation 
for the tidal excitation of modes and discuss, in general, the energy in tidally excited 
modes, the phenomena of tidal resonance, and the effects of rotation and orbital evolution 
on the mode amplitude. In §3 we discuss the relation between the mode amplitude and 
the more observationally accessible quantities such as flux variation, surface velocity, and 
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perturbation to the pulsar's orbit. In §4 we show how the light curve at periastron can 
be used to determine the orbital inclination angle. In §5 we apply our work to two binary 
radio pulsars and to general polytropic stars. Our results are summarized in §6. 



§2 Formalism 



In this section we present our basic method for studying the tidal excitation of modes 
in rotating and non-rotating stars in orbits of arbitrary eccentricity, and discuss the signifi- 
cance of the various modal and orbital parameters. We further present order of magnitude 
estimates for the mode energy. 



Consider a binary system where the primary star, taken to be a point mass, has mass 
M p and its companion star, the secondary, has mass M*. The position of the primary 
star in the orbit is specified by R(t), the separation between the stars, and the azimuthal 



fi* about an axis perpendicular to the orbital plane. 

The perturbation to any physical quantity in the star, for instance the Eulerian density 
perturbation, due to the tidal gravitational field of the primary, can be written as a sum 
over its normal mode amplitudes 



where Sp n £ is the normalized density eigenfunction (normalized such that the energy in 
the mode is unity), Y^ m (#,</>) is the spherical harmonic function, and A n f m is the mode 
amplitude. The pressure (p) and internal gravity (g) modes are uniquely specified by three 
numbers: n, the number of nodes in the radial direction, the spherical harmonic degree 
£, and the azimuthal number m; the surface gravity mode, or the fundamental mode (/), 
has no nodes along the radial direction. We use a shorthand notation, a, to denote the 
collective index (n,£,m). 

The equation for the mode amplitude, A a , can be cast in the form of the following 
second order ordinary differential equation for a forced, damped harmonic oscillator: 



§2.1 Basic equation for tidal excitation &z its general solution 



angle, <f)(t). The secondary is assumed to be rotating as a solid body with angular speed 
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as derived for a non-rotating star in Press and Teukolsky (1977), where 

Qnl= f * drr e+2 8p a (r) (4) 
Jo 

is the overlap integral for the mode, r a is its dissipation rate, R* is the radius of the star, 
and the orbital plane is taken to be at 8 or b = tt/2. It should be noted that our Q n t is equal 
to Press and Teukolsky 's overlap integral divided by u a . This is due to the difference in 
the normalization of our eigenfunctions. Our normalized modes have unit energy, whereas 
Press and Teukolsky's have energy equal to uj 2 a . 

It is useful to define the following reduced forcing function 

f "" {,) = ' <5) 

where a is the semi-major axis of the orbit. The function fi m (t) contains all the time 
dependence of f a (t) and is independent of the mode properties. Since the driving term 
falls off as 1/R e+1 , most of the energy is in the lowest degree modes that can be tidally 
excited, i.e., the quadrupole modes. We rewrite f a (t) in the following compact form: 

fa(t) = fconsKa ftm-, (6) 

where 

f = \ 4n 1 (Ml 

Icons - [( 2 £ + iy-2\ [ Mt 

fi = y/GMt/a 3 is the orbital frequency, and M t is the total mass of the binary system. 
All the modal information is contained in the constant factor K a . For quadrupole modes, 
considered in detail here, f cons is equal to 47r/5 times the fractional mass of the primary 
star. 

Equation (2) can be solved in terms of a Green's function to yield 

A a (t) = f cons ( ) tl 2 f dh exp[-r a (*-*i)] sm[uj' a (t - h)] ft m {h), (8) 

\ U a J J-oo 



K a = L0 2 a Qnl } (7) 



where u' a = u 2 a - V 2 a ^i u a . 

While equation (8) allows us to calculate the mode amplitude for a star with arbitrary 
rotation rate, it does not lend itself to a simple interpretation of the significance of various 
parameters. However, a simple interpretation is possible when the forcing function is 
periodic, such as when the star is non-rotating, and one can analyze the mode amplitude 
using a Fourier series. In the following discussion we assume that the star is non-rotating. 
The effects of rotation on the mode frequency and amplitude are addressed in §2.3. 
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When rotation of the star is neglected, the reduced forcing function fi m (t) is periodic 
with orbital period T, and we decompose it in terms of its Fourier coefficients: 

oo oo 

fim(t) = ]T sin(n0 t) + ]T cos(nfi t). (9) 

n=l n=0 

The Fourier coefficients Cn" 1 ^ and Dn™^ depend only on the eccentricity, e. The mode 
amplitude can easily be calculated using this series and is given below 



A a (t) _ ^ _ sin(nQ t + c/> n ) ^ D n cos(nQ t + c/> n ) 

f cons K a - a ^ [ (n 2_ r 2) 2+4n2c/ 2] 1 /2 A. [ (n 2_ r 2)2 +4n 2 (/ 2] 1 /2 
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where 

tanc/> n = - — — , (11) 

[r z a - n z ) 

and the dimensionless numbers r a and d a are defined below: 

r a = d a = (12) 

It is clear from equation (10) that resonances occur for integral values of r a . Because of 
our choice of normalization for the eigenfunctions (to yield unit energy), the energy in the 
modes, E a , is approximately equal to |A a | 2 . 

For highly eccentric orbits, the tidal force is appreciable only when the star is near 
periastron. Therefore, another relevant time scale is the periastron passage time. We 
define a dimensionless parameter below which expresses the mode frequency in terms of 
this time scale, 

K = ^, (13) 
where Q p is the angular velocity of the star at periastron and is given by 

Up = l (1 _ e) 2 "o- (14), 

The amplitude for a quadrupole mode is equal to (4TT/b)(M p /M t )A a K a (see eq. [10]). 
The dependence of A a on the orbital parameters is entirely contained in A a , which is a 
function of dimensionless parameters, e, r a and d a \ its properties are discussed below. 
In addition, A a depends on the structure of the star, via the parameter K a , which will 
be calculated in §5. The relation of A a to the observed quantities is non-trivial and is 
discussed in 53. 
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The behavior of A a can be understood in terms of the Fourier coefficients. It is easy 
to show that for quadrupole tide -D n 's are real, C n 's are imaginary, and Ci 20 "* = 0. In 
addition, numerical calculations show that lei 22 "*! ~ |I?i 22 ' ) | for large n (n > Q, p /Q, ). 
Moreover, D^°^ falls off as approximately exp(— a e nQ /Q p )/(l — e) 3//2 for n > 5Q p /Q , 
where a e 1.3e -0 ' 25 . The Fourier coefficients for m=±2 are larger compared to m=0 
by a factor of about exp(2a e ). The exponential dependence arises because for eccentric 
orbits, the characteristic time scale for the variation of the tidal force is the periastron 
passage time, O^ -1 ; therefore, the power spectrum for frequencies greater than Q p falls off 
exponentially. 

Tidally excited modes of a star in a circular orbit exhibit oscillations with a period 
equal to that of the orbit and not their natural oscillation periods (so long as the orbital 
evolution timescale is long compared to the mode dissipation time). In the general case 
of an elliptical orbit, the time dependence of the mode amplitude consists of two parts: 
an almost static part that varies on the time scale of the tidal force and a part that 
oscillates with the mode frequency. To be more specific, we can write the reduced mode 
amplitude as A a (t) = A^a\t) + A^^^t), where A^a\t) is equal to the sum of the resonant 
and near resonant terms (n ~ r a ) in the Fourier series and A^a r \t) is the sum of the 
off-resonance terms. The former varies sinusoidally with the mode frequency, while the 
latter is approximately proportional to the tidal forcing function, which peaks strongly near 
periastron. Since the Fourier coefficients fall off exponentially, the oscillation amplitude, 
A^a\t), decreases with r a as exp(— a e r a Q /Q p ) = exp(—a e b a ) for b a > 5. On the other 
hand, because most of the contribution to I« r ' comes from lower order terms in the 
Fourier series, I« r ' is approximately proportional to It should be pointed out that 

this separation of the mode amplitude into two parts only makes physical sense for modes 
with periods comparable to or smaller than the periastron passage time (i.e. b a > 1). 
Modes with b a <C 1 cannot respond quickly enough to the variation of the tidal force near 
periastron. Hence, I« r ' for these modes are not expected to be proportional to the forcing 
function. 

Figure 1 shows the time variation of the reduced mode amplitude (A a ) for three 
different values of b a for the same orbit. For small b a , the resonant term is comparable to 
the sum of the off-resonance terms. Therefore, A a does not vary much from periastron to 
apstron. In contrast, modes with higher b a have large peaks near periastron, arising from 
the static term Aa° r ' ) (t), and much smaller oscillation amplitudes due to the exponential 
fall-off of A^a\t) with b a . The time variation of the reduced mode amplitude for large b a , 
being dominated by the static part, follows closely that of the driving function. Physically, 
this behavior is due to the fact that for large b a , the oscillator adjusts adiabatically to 
the driving function; hence, the mode amplitude is approximately equal to the driving 
function divided by the mode frequency squared. Even the b a = 10 mode does, however, 
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have a non-zero oscillatory component, but its amplitude is a factor of about 2000 smaller 
compared to the static tidal amplitude at periastron. Thus it is exceedingly difficult to 
observe pulsation of modes with large b a values. 

The dependence of A a on r a and d a is simple to understand by making use of equation 
(10). We consider the various possibilities here. When r a is a 'small' integer, that is, when 
r a < 50p/0 or b a < 5, and d a < 1, the r a -th Fourier coefficient is not much smaller than 
the first, and the contribution to the mode amplitude is dominated by the resonant term. 
In this case the oscillatory part of the mode amplitude is much greater than the static 
part, and the mode energy decreases with increasing damping rate as 1/d^. For d a ^> 1, 
several harmonics of the orbital frequency lie within the modal linewidth, and the mode 
energy is proportional to \jd a . 

For 'large' integral values of r a , namely r a > 50p/0 or b a > 5, the Fourier coefficients 
Cii m) and Dii m) are so small (due to the exponential fall-off) that off-resonance terms in 
the series make the dominant contribution to A a . In this regime, the static tide limit, the 
oscillatory part of the amplitude is much smaller than the static amplitude. In addition, 
the mean energy in the tides is essentially independent of the damping rate. 

Finally, for non-integral values of r a , it is convenient to write r a as the sum of the 
nearest integer and Sr a . For d a <C 8r a the mode amplitude is independent of the damping 
rate, and increases as l/8r a . This has the obvious physical interpretation that the ampli- 
tude of an oscillator driven off-resonance is independent of its damping. For d a ^> 8r a the 
behavior is the same as the case of integer r a discussed above. 

~(r) 

The oscillation amplitude is proportional to the product of A a and the overlap inte- 
gral (Q a )- There are thus two opposing effects in the tidal excitation of modes. On the 
one hand, larger n g-modes have smaller frequencies and thus smaller values of b a ; they 

~(r) 

therefore have larger values of A a . On the other hand, Q a are larger for modes with 
smaller n; this condition favors low order modes. For these competing reasons, the modes 
with the largest oscillation amplitude, for most binary systems, are some intermediate or- 
der g-mode (n ~ 5) with 8r a close to zero (see §2.2). It should be noted that the amplitude 
near periastron is dominated by the pi, /, and g\ modes because of their large overlap 
integrals. Since the values of b a are generally large for these modes, the time dependence 
of their amplitudes is the same as the tidal force. Thus we cannot infer much about the 
star from the observation of its light curve near periastron, though it can be utilized to 
determine the orbital inclination angle (see §4). 

To give a feel for the physically relevant range of parameters, we note that the periods 
of low n modes for main sequence stars are of order a few hours to a few days, and the 
damping times are of order 100 years or greater. The energy in these modes due to tidal 
forcing is independent of damping unless their b a value are of order unity and r a happens 
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to be extremely close to an integer (Sr a < d a ). 

§2.2 Order of Magnitude Estimate of Mode Energy 

In this section, we estimate mode frequencies and overlap integrals and use them to 
calculate the energy in tidally excited modes. The frequencies of quadrupolar / and low 
order p-modes are approximately equal to (Christensen-Dalsgaard and Berthomieu 1991) 

uj n2 6 1/4 (n + 1 



where M* and R* are the mass and the radius of the star, and n=0 for the /-mode. The / 
and p-mode frequencies have a weak dependence on the structure of the star which is not 
included in the above expression. Gravity waves can propagate only in convectively stable 
media and thus the (/-mode frequencies are sensitive to stellar structure. For a polytropic 
star of index n p i, and a constant ratio of specific heats, 7, the (/-mode frequencies are given 
below (Christensen-Dalsgaard and Berthomieu 1991) 

fg fGM^ 1 ' 2 

U n 2 



(n + 1) V Rl 

where f g = lb[n p i(^f — 1) — l]/(n p i + 1); the factor of 15 has been obtained by a fit to 
numerically calculated frequencies. 

Since the eigenfunctions are normalized such that there is unit energy in the mode, 
the overlap integral, Q n 2 } is proportional to (Rl/G) 1 / 2 . For n / the wave function is 
oscillatory in the radial direction and this leads to cancellation of the tidal force experienced 
by the mode. Thus the overlap integral decreases with n. We model this dependence 
crudely as (n + and write \Q n 2 \ ~ f q (Rl/G) 1 / 2 /(n + 1)^ . Numerical calculation of 

Q n 2 for polytropic stars with 7 = 5/3 gives f q ~ 1 for / and p-modes, and f q ~ 1/3 
for g-modes. The exponent, /3, for g-modes varies from 1.5 to 3 as the polytropic index 
decreases from 3 to 2 (for an analytical derivation of these results, see Zahn 1970), and 
(3=6 for / and p modes of n p i = l.b stars. 

We estimate A a (t) using equation (10) and find that the energy of oscillations in m=0 
quadrupole modes is 



Ml) [(<Sr„) 2 +<]' 

where n a is the integer nearest to r a . Making use of the above expressions for Q a and 
u a of quadrupole (/-modes and the expression for the Fourier coefficients obtained in the 
previous section, we find the oscillation energy in m = (/-modes to be 



E n2 o ~ E or 



M t 



R± 



d p eri 



f, 



n 



+ 1)2 + 2/9 



exp(-2.66 a e-°- 25 ) 

15 



[(Sr a y +dl](l-e) 
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where d per i = a(l — e) is the periastron distance, and E or b = GM*M p /2a is the orbital 
energy of the binary system. 

Since p-modes have larger frequencies, and thus larger values of b a} than (/-modes, 
p-mode oscillation energies are generally much smaller than (/-mode energies. Moreover, 
we see from equation (15) that the oscillation energy in (/-modes falls off as approximately 
(n + l) -5 for polytropic stars of n p i=3. Thus the mode with the largest oscillation energy 
is typically a (/-mode of moderate n value. The energy in m=±2 modes is larger compared 
to m=0 modes by a factor of approximately exp(5 e -0 ' 25 ). 

If the periastron passage time is large compared to the periods of low order p and 
(/-modes, the tide near periastron is almost static. In this case we can easily estimate the 
mode amplitude by discarding the time derivative terms in equation (2). We find that the 
energy in the tide near periastron, which is dominated by / and low order p and g modes, 
is equal to k E or f,(M p / M*)(R* / d per i) 5 , where k = (2irG/5Rl) Q 2 n2 can be shown to be 
equal to the apsidal motion constant. 

Thus far we have discussed the dependence of the mode amplitude and the average 
energy on the orbital parameters and the mode frequency (through r a , b a and e) as well 
as simple order of magnitude estimates of the dependence on stellar structure. A more 
detailed discussion of the dependence on stellar structure will be given in §5.2. 



§2.3 Effect of stellar rotation on the mode amplitude 



Rotation of the star has been neglected in the previous sections. We now discuss 
the effect of rotation on the mode amplitude. Assuming that the rotation rate is small 
compared to the mode frequency, we can treat its effects using perturbation theory. To 
linear order in the rotation rate, the frequencies of the modes, as seen in the rotating frame 
of the star, are shifted due to the coriolis force and are given by (Unno, et al. 1989) 



m 



m 



'161 



where ^ r nl and £^ are the radial and transverse components of the displacement eigen- 
function (see eq. [34]). Together with equation (8), it is clear that to first order, stellar 
rotation has no effect on modes with m = 0. (For simplicity, we have ignored the per- 
turbation to the eigenfunctions, which could be substantial and would modify the overlap 
integral.) To find the effect of rotation on the amplitudes of tidally excited modes with 
m 7^ 0, we substitute the Fourier series expansion of the reduced forcing function (eq. [9]) 
into equation (8) and carry out the integration which yields the following expression for 
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the reduced mode amplitude 

~ exp(imO t t) 
A n = - 
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It is easy to show that for m = ±2, Re(C 
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tan <£' 



2d a (ms + n) 
~~ ( n + ms) 2 



Im(D 



0, and for n greater than 
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Making use of these in equation (17) we obtain 



A. 



exp(imfi t f 



oo 



D { n m) exp(=Fmfi * T «<MM)) 



J V\j 



n 



m s 



1212 



+ 4C/ 2 . 



'191 



n 



m s 



The mode amplitude observed from an inertial frame is proportional to A a (t) exp(— imQ*t). 
Thus the observed mode amplitude is a periodic function with phase coherence from one 
orbit to another, in spite of the rotation of the star. From the denominator in equation 
(19), we note that resonance occurs when n = r r a + \m\s. Thus, so long as fi* <C uj a} 
the observed mode amplitude is approximately the same as the amplitude of a mode for a 
nonrotating star with frequency equal to u a + fi*(|m| — mJ n t). 

§2.4 Effect of Orbital Evolution on Mode Amplitude Near Resonance 

Orbital evolution causes a resonant mode to move off resonance with time. Thus 
amplitudes of resonant modes do not diverge even in the limit of vanishing damping. We 
showed in §2.1 that the mode amplitude is proportional to l/8r a if d a <C 8r a . Physically 
this is the result of the mode amplitude building up as the mode is kicked in phase at 
periastron over l/8r a orbits. If the change in the orbital period is AT in one orbit, then 
the phase difference of a mode at periastron in two consecutive orbits is uj a AT. Therefore, 
the tidal force will be out of phase with the mode in 2tt / ' ^/u a AT orbits. This implies that 
there is an effective lower limit to 8r a , 8r m i n} for the purpose of calculating the amplitude 
of tidally excited modes. A straightforward calculation yields 



AT 
2tt V~T~ 



I 1/2 



(20) 
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Thus, because of orbital evolution, it is not physically meaningful to consider resonance 
where the ratio of the orbital period to the mode period is any closer to an integer than 
the value Sr m i n . We calculate AT below, which we then substitute into equation (20) to 
determine Sr m i n . 

The decrease in period AT is related to the decrease in orbital energy in one orbit, 
AE or b } as follows 

3 AE or b 



AT 



21 



T 2 E or b 

The change in the orbital energy is equal to the work done on the star by the tidal force, 
which is 



AE, 



orb 



T 



dt / d 3 r v • Fa 



de 



Qnl 



(21 + 1) 



T 



dt 



Y e * m (Tr/2,<f> orb )dA n 



£m 



dt 



(22) 

where v is the velocity induced in the star by the tidal force and is equal to ^ a £ a dA a /dt. 
We calculate the time integral in the above equation using equations (9) and (10) and the 
result, for b a ^> 1 for the low order modes, i.e. the static tide limit, is given below 



AE : 



orb 



E, 



orb 



32nk h(e) 



de 



M p M t 
Ml 



R* 
a 



(23) 



where k is the apsidal motion constant given by 

2ttG \ - 



k 



5RI ^ 



Ql 



25 



(24) 



Ttide is the dissipation rate of the tide which is in general different from the mode dissipa- 
tion rate introduced earlier and 



Me) 



3072 + 47616e 2 + 97920e 4 + 35520e 6 + 1200e 8 



2048(1 



!)15/5 



(25) 



In deriving equation (23) we took T q /to^ = T tide(R% I GM*) . This is a crude description for 
the dissipation of static tides in terms of the normal modes of the star and corresponds to 
assuming that the lag-angle of the star associated with different modes is equal. (It is easy 
to show that the lag-angle associated with a mode, for slowly rotating stars, is 2r g Oo/^a)- 
We note that our expression for orbital energy change (eqs. [23]-[25]) is identical to that 
of Hut (1981) for nonrotating stars. 

Combining equations (20), (21) and (23) we find the minimum possible value of 8r a 
to be 



(fr rara ) 2 =24r a k h(e) 



de 



M p M t 
Ml 



R* 
a 



(26) 
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The lag-angle of the tidal bulge, 8 t ide, is related to T t ide- $tide ~ ^tide^o{Rl/ GM*). This 
can be used to express the above limit in terms of the lag-angle. 

The angular momentum of the orbit also changes due to tidal interaction, which causes 
a change in the orbital eccentricity. However, in the absence of energy transfer from the 
orbit to the star the orbital period remains unchanged. Therefore, a change in e alone 
does not change the phase relationship of a mode and the forcing function at periastron. 
Thus we do not expect our result given in equation (26) to be significantly modified when 
angular momentum transfer is included. Finally, we note that it can be easily shown that 
the apsidal motion of the system has the same effect on the mode amplitude and resonance 
as the stellar rotation discussed in §2.3. 

§3 Observation of pulsation 

The pulsation of a star can be detected by one of the following three techniques: (1) 
photometric or flux variations; (2) spectroscopic observation, i.e. , by measuring the time 
dependent surface velocity; or (3) for binary pulsar systems, the pulse arrival delay due to 
the periodic variation of the neutron star orbit caused by the time dependent quadrupole 
moment of the oscillating star. The mode amplitude calculated in the last section needs to 
be transformed to yield these observational quantities. The details of these transformations 
are described below. It should be noted that one can distinguish tidally induced pulsation 
from the intrinsic pulsation of the star because the former has a definite phase relation to 
the orbital motion. 

§3.1 Flux variation associated with pulsation 

Associated with the tidally excited modes there are variations of the star's luminosity. 
The observed flux variation arises from the change in the surface temperature of the star 
as well as the change in its projected surface area. To first order, we can calculate these 
two effects separately and add them to obtain the total flux variation. 

The flux change at the stellar surface, associated with a mode (n,^,m), is given by: 
8F ntm {R*,6, 4>) = A ntrn 8F nt {R*)Y trn {6, </>), (27) 

where 8F n t(r) is the flux variation eigenfunction. To calculate the observed flux variation, 
^rUmi we nee d to integrate the above expression over the hemisphere of the star facing 
us. Thus the expression for the fractional observed flux variation is 

^nlm = J fi f LD {fi)Y im {e,<f ) )8F ni {R,) 

F° bs fdSliif LD (ii)F ' 1 ) 

where 

/j = n s ■ n obs , (29) 
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n s is a unit vector normal to the stellar surface, n & s is a unit vector pointing to the 
observer (from the center of the star), F is the unperturbed flux, and f LD is the limb 
darkening function which can be written as 



Ild (aO = a o + aiji + a 2 fi' 



(SO) 



where a 0.3726, a\ ~ 0.6500, and a 2 ~ —0.0226 (see, e.g., Kopal 1959). The expression 
in equation (28) does not include the variation in the projected area of the star, which is 
calculated below. 

Specializing to i = 2 modes and making use of the last four equations, we find the 
observed flux for the m=0 modes to be 



3F°2Q _ fconsO LD / 5 

pobs ~ 4 V 167r 



Kn2$F n 2(R* 

F 



A n 2o [3 



cos 2 i 



(SI) 



and the observed flux for the sum of the m=— 2 and m=2 modes to be 



$Fn22 _ fconsa LD / 15 

F ohs ~ 2 V 32tt 



Kn2$F n 2(R* 

F 



Re(A n 22) cos(2c/)o — 20*^ 
Im(A n 22) sin(2(/> — 20*^ 



f32l 



sin 2 i. 



where i is the inclination angle of the orbit, and <f>o is the azimuthal angle of the line-of- 
sight projected onto the orbital plane and is related to the longitude of periastron (u) by 
(f>o = io — 7r/2. The factor a LD , which arises from limb darkening, is defined as 



15a + 16ai + 15a2 
5(3a + 2a\ + 1.5a2) 



(331 



Note that the observed flux variations for m=0 and m=±2 modes (eqs. [31] & [32]) have 
a different dependence on the inclination angle. This can be used to determine the orbital 
inclination angle (see §4 for details). 

The calculation of the flux variation eigenfunction, SF n £ } requires solving coupled 
linear differential equations for nonadiabatic oscillations, which is quite involved. However, 
for most purposes it is reasonably accurate to use the relation SF n 2(R*) / F ~ 4AT(r = 
2/3)/T, where AT(r = 2/3) is the Lagrangian temperature perturbation eigenfunction 
at optical depth r=2/3 which can be calculated using an adiabatic oscillation code. The 
product, K a 8F n t(R*), which contains the dependence of the observed flux variation on 
the structure of the secondary star, will be calculated for several different stars in §5. 

We next calculate the change in the projected surface area of the distorted star as- 
sociated with different modes. The change in the area also causes the observed flux to 
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vary. The distortion of the surface due to one mode can be found from the displacement 
eigenfunction, which can be written as 

where and £^ are the radial and transverse displacement eigenfunctions, respectively. 
These eigenfunctions are normalized so that the energy in the mode is unity, i.e., 



(35) 



The fractional change in the projected area of the star associated with a mode, 
^S n £ m /irR^, can be expressed in terms of the surface displacement. The fulx variation 
arising from this change in projected area can be shown to be identical to the expression 
for 8F° h 2 s m lF ohs of equat ions (31) and (32), provided that we replace a LD [SF n 2 {R* ) / F] in 
these equations by 

' 8& 2 (fl»)-6&(fl») ' _ 

This allows us to define an effective flux variation eigenfunction SF^[^ as 

"4£ 2 (i2„)-3&(i2„ 



F F 



R± 



(36) 



The total observed flux variation associated with modes is obtained from equations (31) 
and (32) by replacing 8F n 2/F in these equations with the effective flux eigenfunction as 
given by the above equation. 

§3.2 Surface velocity amplitude associated with modes 

The calculation of the surface velocity proceeds in a similar manner to the flux vari- 
ation calculation. The line-of-sight surface velocity, integrated over the stellar hemisphere 
facing the observer, is given by 

Vntm = j dtl \n s ■ n obs \ n obs • * n/ro (i2„), (37) 

where the first factor in the above integrand weighs the velocity by the projected area. 

Carrying out the integral given in equation (37) for £ = 2, we obtain the following 
expression for the observed velocity amplitude for an m=0 mode 



Vn°20 = fcons\\j^U n2 K n2 [C n2 (R*) + 3&(i2,)] A n20 [3 cos 2 i - 1] (38) 
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and for the sum of the m = ±2 modes 



Trobs f 

V n22 — Jco 



nsM , K -w n2 ir n2 [Cn2(R*) + 3£* 2 (#*)] sin 2 i Re(A n22 ) cos(2</> - 2fi*t) 

107T L 



— Im(A n22 ) s'm(2(f> — 2flJ) 

J (39) 

These expressions have the same functional dependence as the flux variation derived in 
the previous subsection. 



§3.3 Perturbation to pulsar orbit due to pulsation 

The contribution to the gravitational potential at the position of the primary star due 
to one mode of pulsation of the secondary can be easily calculated using equation (1) and 
is given below 

4TrGQ a Y em (Tr/2,<f> orb (t))A a (t) 



8^, 



(40) 



(2£ + l)i^+! 

The equation for the perturbation to the primary's position, Ri, due to this potential, to 
first order in Ri , is as follows 



dt 2 



[Ro/af 



Ri — 3Rq(Rq ■ R-i 



(41) 



where Ro is the primary's unperturbed position with respect to the center of mass, Ro = 
Ro/|Ro|, and M t is the total mass in the binary system. The first term in the above 
equation arises due to the change in the gravitational force on the primary star because of 
the perturbation to its position. In general the solution of equation (41) must be obtained 
numerically. However, for modes of small b a the mode amplitude A a is a sinusoidal 
function, and if r a ^> 1 (highly eccentric orbit), the first term on the right side of the 
above equation can be ignored, and equation (41) can be solved to yield the ratio of the 
perturbed to the unperturbed velocity of the primary. The result below is for quadrupole 
modes 



lu^Ril 167T 2 { {l LO a Q 2 a M p 



a0 



25 



M t M, 



A a (t)Y em (Tr/2,<f> orb (t)) 
a 2 (i?/a) 4 



(42) 



Using results from §2.2, we find that the fractional perturbation to the primary's orbital 
velocity is typically very small. 

It is useful to express the perturbation to the orbit as a delay in the pulse arrival time, 
8t a , when the primary is a pulsar. This delay is a periodic function and is given below 



8t c 



sin i (Ri • n )/ c, 
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where n , a unit vector, is the projection of the observer's direction on the orbital plane. 



We describe a technique for determining the inclination angle of an eccentric orbit 
by fitting the observed light curve near periastron to theoretical light curves. The main 
physical idea behind this method is simple. In the static tide limit, when periods of low 
order modes are small compared to the periastron passage time, the star adjusts its shape 
to the instantaneous tidal force and thus the tidal bulge points in the direction of the other 
star irrespective of its rotation rate. Therefore, as the secondary star makes its passage 
through periastron, the tidal bulge continues to point toward the primary star, and the 
time dependence of the tidally distorted surface of the star projected onto the plane of the 
sky gives rise to the observed light curve. Since the observed part of the surface of the star 
depends on the inclination angle, the shape of the light curve also depends on the angle i. 

For an orbital inclination angle, i } and a longitude of periastron, u, the observed flux 
variation, SF b S} is the sum of the contribution from m = ±2 and m = modes and can 
be written in the following form (see eqs. [31] and [32]) 



8F hs{t) oc sin 2 i [sin(2cu — 2fi*t)\E , 2 (t) — cos (2w-2ft**)#i(*)] + (3cos 2 z - 1)# (<), (43) 



where \I/ 1 and are the real and imaginary parts of the mode amplitude summed over all 
m = 2 modes (weighted by K^fiF^^), as seen in the rotating frame of the star, and ^>o 
is the sum over all m = modes. In the limit of static tide, the tidal bulge always points 
towards the other star; therefore, the observed flux variation is independent of the stellar 
rotation rate fi* (see the derivation below). We see from equation (43) that the shape of 
the observed light curve in general depends on the orbital inclination angle because of the 
different relative contributions from m = and m = ±2 modes. For the particular case 
of circular orbits, however, the amplitudes of m=0 modes are time independent so that 
the shape of light curve is independent of i. Of course, the magnitude of the observed flux 
variation still depends on i } but because of the uncertainty in stellar mass and structure 
this is of limited use in determining i. 

If the dimensionless mode frequency, b a , for the first few p and (/-modes is greater than 
about 5, then the time derivative terms in equation (2) can be neglected and so the mode 
amplitudes are equal to f a /(jj^. In this static tide limit the function ^o(t) 1S proportional 
to the reduced forcing function /20, and ^i(t) & ^2(t) are proportional to the real and 
imaginary parts of ^22 (it is easy to show that the tidal bulge points in the direction of the 
other star in this case). Thus using equations (5) and (43) we obtain the following explicit 
expression for the light curve 



§4 Determination of orbital inclination angle 



SF obs (t) 



F c 



1 — 3 sin 2 i sin 2 ((f> rb (t) 




(44) 



(R(t)/a) 
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where F c is a constant factor that depends on the stellar mass and structure and is inde- 
pendent of the orbital parameters. F c can be calculated either by determining the normal 
mode eigenfunctions and the associated flux variation at the surface or using von Zeipel's 
theorem (cf. Kopal 1959), according to which the emergent radiation flux, in the static 
tidal case, is proportional to the local gravity. We note that the rotation rate of the star 
has dropped out from the above equation for the observed flux variation as was claimed 
earlier. Thus a comparison of the shape of the observed light curve with the function given 
in equation (44) provides a robust method, which does not depend on the unknown mass, 
structure and rotation rate of the star, for determining the orbital inclination angle. 

Light curves, corresponding to several different inclination angles, for a 10 Mq main 
sequence star in an orbit of e=0.4 and period of 100 days are shown figure 2. The value of </>o 
is chosen to be 0°. We emphasize that the magnitude of the flux variation in these graphs 
depends on the mass and structure of the star as well as the inclination angle. However the 
shape, which is what we suggest should be used to compare with the observed light curve 
to determine the inclination angle, does not. We note that for small i } the contribution 
to the observed light curve is dominated by m=0 modes. As i increases, the contribution 
from the m = ±2 modes increases, resulting in a significant change in the shape of the 
light curve. The shape of the light curve is most sensitive at intermediate values of i 
(20 < i < 70); however, there are quite noticeable variations for all i. 

The order of magnitude of the flux variation can be estimated crudely from the change 
in the Lagrangian temperature at the point in the star's atmosphere where the optical depth 
is 2/3. Since the tide is very nearly static, we can integrate the hydrostatic equilibrium 
equation from the surface over a fixed column of gas corresponding to an optical depth of 
2/3, with and without the tidal force, to obtain the fractional change in the Lagrangian 
pressure and find it to be Sg/g + F t /g = 28R*/R* + 2(M p / M*)(R* / d) 3 , where g is the 
gravitational acceleration at the stellar surface, F t is the tidal acceleration, d is the distance 
between the stars, and 8R* is the change in star's radius, which is approximately equal to 
(M p /M*)(R* I d) 3 R*. Thus the fractional change in the luminosity of the SMC star, for 
which Mp/M* ~ 1/7 and R*/d ~ 0.2 at periastron, is estimated to be about 0.5% or 5 
m-mag. Our numerical value is a factor of about 2 larger (see §5.1). 

The so called ellipsoidal variable stars are believed to undergo flux variations due to 
tidal distortion which are correlated with the orbital phase (cf. Kopal 1959, Avni & Bah- 
call 1975, Morris 1985). It is well known that the magnitude of the observed light variation 
of the ellipsoidal variable stars depends on the orbital inclination angle. However, it is not 
appreciated that the shape of their light curve also depends on the orbital inclination angle 
if the orbit is non-circular. Finally, we note that our proposed method for determining 
orbital inclination angle is expected to work best in those systems where one of the stars 
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is compact so that the interpretation of light curve is not complicated by eclipses. 

§5 Applications 

In this section, we apply our formalism to two radio binary pulsar systems, PSR 
J0045-7319 and PSR B1259-63, and discuss the possibility of direct observation of the tidal 
excitation of modes in these systems. We also present results of numerical calculations for 
general polytropic stars in order to further examine the dependence of tidal excitation on 
stellar structure (§5.2). 

§5.1 Applications to PSR J0045-7319 and PSR B1259-63 

Recently, Kaspi et al. (1994) have reported the discovery of a radio pulsar (PSR 
J0045-7319) in the Small Magellanic Cloud. Measurements of the Doppler time delay have 
indicated that the pulsar is in a highly eccentric (e = 0.80798) orbit with a period of 51.169 
days. Kaspi et al. (1994) have optically identified the companion star as a 16th magnitude 
B-star with mass of lOM©. In our calculations, we also assume that the mass of the 
pulsar is 1.4M0. 

Using a B-star model kindly provided to us by the Yale group, we have computed 
the eigenfrequencies and the overlap integrals (Q n i) for quadrupole oscillations. We find 
that the /-mode frequency for this system corresponds to a b a value of approximately 
21.7, and the g$ mode has a b a value of 3.96. The frequencies of m = ±2 modes are 
shifted by approximately ±l/iHz, if the B-star is rotating synchronously at periastron, 
which corresponds to a shift in the value of r a of about ±4.4. The radiative damping time 
depends on the mode order: for g2 and (/3-modes we estimate it to be about 10 3 years, and 
for (/10 mode to be about one year. The turbulent damping time for these modes is a factor 
of a few larger. The radiative damping time for the P2 mode is about a year. Thus for 
none of these modes do we expect the mode energy to have any dependence on damping. 

For the lOM© ZAMS stellar model the product K n 2 SF n 2(R*) / F is very nearly inde- 
pendent of n for (/-modes of order 4 to 10, beyond which this product starts to decrease. The 
probability that one of the approximately 25 low to intermediate order g-modes (counting 
m=0, and ±2 as three distinct modes due to their different frequencies because of rota- 
tional splitting) has a value of r a that is close to resonance (|^r a | < 0.05) is very high. 
In fact, for the stellar model we use, the g-j mode has 8r a of 0.02, which is the smallest 
value for the first 10 harmonic g modes. The smallest physically meaningful value of |$r a |, 
resulting from the orbital evolution considerations (see §2.4), is estimated using equation 
(26) to be 0.003. Thus, Sr a =0.02 for the g-j mode is physically acceptable. The g-j mode, 
which has the largest amplitude, has a period of 1.07 days (corresponding to b a = 2.96). 
The stellar flux variation associated with this mode is 2.3 m-mag, and the surface velocity 
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amplitude is 70 m s 1 (calculated using eqs. [38] and [39]). For comparison, the flux 
variation associated with the /-mode, away from periastron, is about 10 //-mag. 

The perturbation to the pulsar orbital velocity due to the time dependent quadrupole 
moment of the g-j mode has been calculated by integrating equation (41) and is found to 
be 10 -3 cm s _1 . The density fluctuation in the convection zone of the star also causes a 
random perturbation to the pulsar's orbital velocity that is estimated to be about 10 -6 



It should be stressed that the amplitudes of these tidally excited modes (shown in 
Fig.l) have perfect phase coherence from one orbit to another which should be helpful in 
observing these small amplitude oscillations. On the other hand, since the mode with the 
largest observed amplitude is one with a frequency closest to resonance, it would be very 
difficult to make a precise identification of the order (n) of the observed mode, rendering 
its potential use in determining the properties of the star rather limited. 

The light curve for PSR J0045-7319 near periastron, dominated by / and low order p- 
modes, is shown in figure 3 for several different orbital inclination angles. The magnitude 
of the observed flux variation is about 10 m-mag and is calculated assuming that the 
secondary star has a mass of 10 M©. The uncertainty in our estimate of the flux variation 
is about 20%, reflecting an uncertainty in the mass of the B-star of 2Mq. The shape of 
the light curves in figure 3, however, is independent of the mass of the star, and can be 
used to determine the inclination angle of the orbit. 

We have calculated the periastron advance of the PSR J0045-7319 system, due to 
dynamical tides, by integrating equation (36) and summing the result over all modes. The 
advance is found to be 0.7 arcsec per orbit. This result is in excellent agreement with the 
classical apsidal motion formula (Claret & Gimenez, 1993); our result is about 1% larger 
compared to the classical formula. Thus the static tide approximation for calculating 
apsidal motion is an excellent approximation for this system. 

The recent accurate observations of PSR B 1259-63 (Johnston et al. 1994) indicate that 
its eccentricity and orbital period are approximately 0.8698 and 1236.79 days, respectively. 
The companion is a Be-star with a mass of 10M©. Thus the periastron distance is 
~ 20i?*, and the b a values for the modes of the Be-star are about a factor of 13 larger 
than those for the SMC B-star. Therefore, the amplitude near periastron, which falls off 
as the distance cubed, is expected to be a factor of about 70 smaller than that of the SMC 
system, and the oscillation amplitudes of even the g modes, because of the large b a values, 
are too small to have any observational consequences. 



In this section, we apply our discussion of mode amplitude to a general polytropic 




§5.2 Application to Polytropic Stars 



20 



star. In particular, we examine how the flux variation and surface velocity depend on the 
structure of the star. Our calculations include the /-mode as well as low order p-modes 
and (/-modes for polytropic stars with index ranging from n p i = 1.5 to 3. 

For a polytropic star, the eigenfrequencies and overlap integrals can be expressed in 
nondimensionalized forms which depend only on the index n p {. We have calculated these 
quantities without the Cowling approximation, and our results are in excellent agreement 
with values quoted in Robe (1968) and Lee and Ostriker (1986). 

The dependence of the mode energy on the stellar structure is contained in the pa- 
rameter K\. We calculate this numerically as a function of polytropic index, n p i, for the 
quadrupole /-mode and low order p-modes and (/-modes. The results are plotted in figure 
4 in the units of GM^/R*. Mode energy is obtained by multiplying K 2 2 , \A n 2 m \ 2 (cal- 
culated using eq. [10]), and ff. ons - While the /-mode and the p-modes have comparable 
values of K 2 , the (/-modes are few orders of magnitude smaller. This is a result of two 
separate effects. First, the overlap integral is largest for the /-mode and decreases with 
increasing mode order n. Second, the eigenfrequencies for the p-modes are much larger 
than those for the (/-modes. Thus, even though p-modes and (/-modes could have compa- 
rable overlap integrals, K a still differs significantly. Figure 4 shows that K\ increases with 
n p i for all modes. This might suggest that higher n p i polytropes are more susceptible to 
tidal excitations. However, it is important to keep in mind that the mode energy is also 
proportional to |A a | 2 , which has a sensitive dependence on the mode frequency through 
b a ; \A a \ 2 decreases exponentially with b a . Since for a fixed mode, the mode frequency (and 
hence b a ) tends to increase with n p i, A a tends to decrease with n p {. Therefore, depending 
on the orbital parameters, the actual mode energy might in fact decrease with n p {. 

The fractional flux variation from the star depends on the dimensionless parameter 
K a 8F n t(R*) I 'F ', which is approximately equal to 4A' Q ,AT(r = 2/3)/T, where T(r = 2/3) 
is the Lagrangian temperature variation at optical depth of 2/3. We calculate AT numer- 
ically using an adiabatic oscillation code. This is then used to compute K a 8F n t(R*) / 'F ', 
and the results are shown in figure 5 for several different quadrupole modes as a function 
of n p i for a fixed stellar mass of 1.0 M©. The p-modes correspond to acoustical compres- 
sion or expansion, while the (/-modes are effectively incompressible. Hence, p-modes have 
correspondingly larger values of AT/T. For polytropic stars, the dimensionless parameter 
K a 8F n t(R*) I F is found to be almost independent of the mass of the star. 

The surface velocity amplitude depends on the parameters K a u a il r a and K a u a i^ (see 
eqs. [38] and [39]). These parameters are shown in figure 6 in the units of ^GM^/R*. 
The p-modes have predominantly radial displacement at the surface, while (/-modes have 
predominantly transverse displacement. 

The shift of the mode frequencies due to rotation can be calculated using the dimen- 
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sionless number J n 2 (see eq.[16]) which is shown in figure 7, for several different modes, as 
a function of n p i. For p-modes, J n 2 approaches as n increases since the transverse dis- 
placement becomes negligibly small compared with the radial displacement. For (/-modes, 
it is the radial displacement that becomes negligible as n increases; hence J n 2 approaches 
l/£{£ + 1) or 0.1667 for quadrupole modes (see, e.g., Unno et al. 1989). Therefore, J n 2 is 
bounded above by the /-mode and is always less than 0.5. 

In order to obtain the actual flux variation or surface velocity, one needs to multiply 
the relevant products (displayed in Figs. 5 and 6) by the reduced mode amplitude A a , 
which depends on the orbital and modal parameters through r a and e. As an example, 
consider a binary system with e = 0.7, consisting of a 1.0 Mq compact primary, and a 5.0 
Mq secondary star of polytropic index 2. To find the flux variations associated with the 
pi and (/3-modes of the secondary, we read off from figure 5 that the values of K a 8F n tjF 
are approximately 100 and 0.1 respectively. Assume that the orbital period is such that 
the gs mode has r a =23.8 (6 a =3), and the p\ mode has r a = 188.5 (6 a =23.8). Using the 
procedures outlined in §2, we calculate A a for the gs (pi) mode with m=0 and find the 
peak value of the reduced amplitude at periastron to be 0.02 (3.3 X 10 -4 ) and an oscillation 
amplitude of 0.01 (1.3 X 10 -7 ) away from periastron. The product of A a , K a SF n t >/ F , and 
fcons (=2-1) yields the flux variations for these modes to be 20 (33) milli-magnitudes at 
periastron and 10 (0.013) milli-magnitudes for the oscillations away from periastron. 

For concreteness, we have calculated in figure 8 the fractional luminosity variations, 
SF°22/F, for several low-order (/-modes in a n p i = 3 polytrope as a function of orbital 
period. To prevent close resonances from confusing the general trend, we have modified 
the mode frequencies slightly such that 8r a = 0.1 for all the modes considered. The mass 
of the secondary polytropic star is taken to be 5 Mq, and the mass of the primary star is 
2 Mq. To illustrate the dependence on eccentricity, e has been chosen to be 0.8 and 0.4 
(Fig. 8a and Fig. 8b, respectively). Figure 8 shows that SF^^/F falls off rapidly with the 
orbital period and that it is significantly larger for more eccentric orbits. For larger orbital 
periods, variations due to higher overtone (/-modes dominate because of their smaller b a 
values. 

§6. Summary and discussion 

We have calculated the energy and amplitude of the tidally excited modes of a star 
in a binary system of arbitrary eccentricity. The effects of damping, resonance, orbital 
inclination, and rotation of the star on the observed mode amplitude have been investigated 
in detail. The important parameters for this problem are: b a , the ratio of the periastron 
passage time to the mode period; 8r a , the ratio of the difference between mode frequency 
and the nearest resonance frequency (resonance frequencies are integer multiples of the 
orbital frequency) to the orbital frequency; and d a , the ratio of the orbital period to the 
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mode damping time. 

Modes of large b a couple poorly to the tidal force. The oscillation amplitude decreases 
with increasing b a as approximately exp( — 1.3&Q,e -0 ' 25 ), where e is the orbital eccentricity. 
The p, /, and low order (/-modes of main sequence stars, for most binary systems, have 
values of b a that are much greater than 1; thus, their oscillation amplitudes are exponen- 
tially suppressed. Since the overlap integrals decrease with increasing mode order (n), the 
moderate order quadrupole (/-modes with b a ~ 1 and the smallest value of 8r a , i.e., the 
modes most nearly resonant, have the largest oscillation amplitudes. 

The mode amplitude is independent of damping so long as d a is less than 8r a . This 
condition is expected to be satisfied for the low to moderate order g modes of main sequence 
stars which have damping times of at least a few years. In the limit of d a <C 8r a the mode 
amplitude is proportional to l/8r a , as expected of an oscillator driven off-resonance. Due 
to orbital evolution, a mode which is in perfect resonance with the orbit, i.e. 8r a =0 } moves 
off resonance with time. This places a severe upper limit to the amplitudes of modes near 
resonance (§2.4). 

The mode amplitude and phase of tidally excited modes is the same from one orbit 
to another. This property should be useful in observational searches for tidally excited 
modes as it distinguishes them from intrinsically excited modes. Rotation of the star (in 
the limit of slow rotation rate) causes the frequencies of m ^ modes to be shifted, but 
otherwise does not change the properties or time dependence of the mode amplitude. 

We have applied this tidal excitation theory to the recently discovered pulsar in the 
SMC, PSR J0045-7319, which is in a binary system of eccentricity 0.81 and orbital period 
of 51 days. The companion star, believed to be a 10M© B-star, has a periastron distance of 
about 5 stellar radii. Thus this system is a good candidate for studying the tidal excitation 
of modes. Using a 10M© main sequence stellar model, we find that the mode with the 
largest oscillation amplitude is the g-j mode, which has 6 a =2.96 and 8r a =0.02. We find 
that the flux variation associated with this mode is 2.3 milli-magnitude, and the surface 
velocity amplitude is 70 meter s _1 . The perturbation to the orbital velocity of the pulsar, 
as a result of the time dependent quadrupole moment of the B-star due to this mode, is 
about 10 -3 cm s _1 at periastron. These amplitudes are independent of mode dissipation 
time because the value of d a is much less than 8r a (the radiative and turbulent dissipation 
time scales are about 100 years). 

It should be pointed out that the small value of 8r a for the g-j mode is a coincidence. 
A change of 0.2% in the mode frequency changes 8r a by 0.1. Thus 8r a should be regarded 
as a random variable. Therefore, given 25 modes that are similar in terms of their coupling 
to the tidal force, we expect one of them to have \8r a \ less than 0.02. In fact, modes g^—gio 
of the Yale star model have very similar mode amplitudes for a given 8r a . Thus it is not 
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surprising to find one of them, counting m=0, ±2 as three different modes because rotation 
lifts the m-degeneracy of the frequencies, with 8r a of 0.02. Orbital evolution places a lower 
limit on |$r a |, which we estimate to be 0.003. 

The tidal distortion of the star near periastron is dominated by low order quadrupole 
p and (/-modes and the /-mode, for which b a is much greater than one, and therefore the 
light curve near periastron follows the time dependence of the tidal force closely. Since 
the observed light curve depends on the projection of the tidally distorted star onto the 
plane of the sky, which in turn depends on the inclination of the orbit, the shape of the 
light curve near periastron allows the determination of the inclination angle of the orbit for 
eccentric orbits. In particular, this technique should be applicable to the SMC pulsar, PSR 
J0045-7319, for which we estimate the flux variation near periastron to be 10 m-mag, 
and the surface velocity to be about 30 m s _1 . 

For the SMC pulsar, we further determine that the periastron advance due to the 
tidal distortion of the B-star, treated dynamically, is about 0.7 arcsec per orbit. This is 
in excellent agreement (ss 1.0 %) with the static tidal calculation of the classical apsidal 
motion formula. The precession due to general relativity is estimated to be 2.0 arcsec per 
orbit and that due to the quadrupole moment due to rotational distortion (for synchronous 
rotation) is about 2.2 arcsec/orbit. Thus the total apsidal motion compares reasonably 
well with the current observational value of about 5 arcsec per orbit (Kaspi et al. 1994). 

For the only other known binary radio pulsar with a main sequence companion, PSR 
B 1259-63, the flux variation near periastron is approximately 70 times smaller than that 
of the SMC system, and the oscillation amplitude of all g and p-modes is negligibly small. 

We have also calculated the stellar structure dependence of the mode energy, flux vari- 
ation, and surface velocity amplitude by considering polytropic stars of index n p i ranging 
from 1.5 to 3. These quantities can be written as a product of the two terms: a term that 
only depends on the reduced mode amplitude (A a ) and a term that depends on the stel- 
lar structure (K a ,8F a / F, etc.). The stellar structure dependent term generally increases 
with n p i. However, this does not necessarily mean that higher n p i stars give rise to larger 
tidal effects. The mode frequency (and thus b a ) increases with n p i, and consequently A a 
decreases with n p i. Therefore, depending on the orbital parameters, tidal effects might in 
fact be larger for lower n p i stars. 
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Appendix 



In this appendix, we estimate the perturbation to the primary star's orbital velocity 
due to density fluctuations in the convection zone of the secondary star. 

We first calculate the quadrupole moment of the star due to fluctuating eddies in 
one scale height thick layer (H) of fluid at a radial distance r from the center. Let the 
mean mass of this layer be Mh and the magnitude of the fluctuation be 8Mh- Since the 
fractional density fluctuation, Sp/p, associated with turbulent flow is M 2 H , where Mh 
is the Mach number of the turbulence, therefore the fractional mass fluctuation of this 
layer is 8M H /M H ~ HM 2 H /r. Let us assume that the luminosity carried by convec- 
tion is /X*, where L* is the luminosity of the star. According to mixing length the- 
ory fL* ~ 4irpr 2 C 3 M 3 H ~ MhC 3 M 3 h /H, where C ~ \/gH is the sound speed and 
g is gravitational acceleration. Thus the RMS fluctuation of the quadrupole moment 
8Q H ~ r 2 8M H ~ rM^ifL^Hfl 3 /g. Since 8Q H oc (i^Mjj) 1 / 3 oc (pH 3 ) } the contribu- 
tion to the quadrupole moment is dominated by eddies at the bottom of the convection 
zone. So from here on we will take H to be the scale height at the bottom of the convection 
zone. The perturbation to the gravitational potential due to the quadrupole moment 8Qh 
is 



5V 



4ttG8Q h 

5R 3 



We can determine the perturbation to the primary's orbital velocity by solving equation 
(41) with Sty given in the above equation. If the convective time scale is much smaller 
than the orbital period, so that orbital modes are not resonantly excited by convection, 
we can discard the first term in equation (41), and the solution of the resultant equation 
is given below 



a 



c/Ri 



dt 



~R 



a 
R 



tint 



o^kh 



where f r =r/R* } tn=H /Vh is the convective time, Vh is the convective velocity, and 
tkh=GM 2 / (L*R*) is the Kelvin-Helmholtz time scale of the secondary star. The dimen- 
sionless factors / and f r are roughly of order 1. Therefore the random fluctuation to the 
primary's orbital velocity, due to a lOM© main sequence secondary star in an orbit with 
a/R* ~ 25, is a factor of 10 13 smaller than the mean orbital velocity. 
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Figure Captions 



FIG. 1. — The reduced mode amplitude A a as a function of time in one orbit, shown for 
several values of b a for m=0. Here we take e = 0.8, 8r a = 0.1, and d a = 0.01. t = 
corresponds to periastron. Note that the oscillation amplitude falls off exponentially with 
b a , whereas the peak at periastron falls off more slowly. For b a = 2, the oscillation ampli- 
tude is comparable to the periastron amplitude. For b a = 10, the oscillation amplitude is 
so small that the oscillations can hardly be seen in Fig. lc. In Fig. Id, a portion of Fig. 
lc near apstron is expanded in scale to show the oscillations more clearly. These tidally 
excited oscillations have phase coherence with the orbital motion; as such, they should 
be readily distinguished from intrinsic oscillations of the star. The m = ±2 modes have 
similar behaviors but tend to have larger oscillation amplitudes for a fixed b a . 

FIG. 2. — Light curve near periastron for a 10 Mq star in an orbit with e = 0.4 and T = 
100 days for several values of the orbital inclination angles i; the unit of flux variation is 
arbitrary, and t = corresponds to periastron. The value of </>o is chosen to be 0°. The 
shape of the light curve varies as a function of i due to the relative contributions from the 
m = and m = ±2 modes, which project differently onto the plane of the sky. While 
the magnitude of these curves depends on the mass and detailed structure of the star, the 
shape does not. Hence it can be used to determine i even when the stellar structure of the 
star is poorly known. 

FIG. 3. — Light curve near periastron for the B-star in the SMC pulsar system (PSR J0045- 
7319) for several orbital inclination angles. </>o is taken to be 25.24°, corresponding to a 
longitude of periastron of u = 115.24°, as quoted in Kaspi et al. (1994). The shape of 
these curves can be used to determine the inclination angle for this system. 

FIG. 4. — K 2 a in the units of GM^/R* as a function of polytropic index, n p i, for the 
quadrupole /-mode and low order p and (/-modes. The parameter K\ contains the depen- 
dence of the mode energy on stellar structure. Note that K\ for (/-modes is a few orders of 
magnitude smaller than for p-modes even though their overlap integrals are comparable. 
This is because the eigenfrequencies for the (/-modes are much smaller than those for the 
p-modes. 

FIG. 5. — \K a 8F n fJ F\ as a function of polytropic index, n p i, with stellar mass fixed at 
IMq. This dimensionless parameter contains the dependence of the fractional luminosity 
change on the stellar structure. Since p-modes are compressible whereas (/-modes are 
effectively incompressible, the p-modes have larger values of the Lagrangian temperature 
perturbation, AT/T, leading to larger values of K a 8F n tjF. Furthermore, this parameter 
varies monotonically with n p i for p-modes, but has a minimum for (/-modes corresponding 
to a node at the surface for AT. 
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FIG. 6. — \K a Lo a ^ r a \ and \K a uj a ^\ in the units of \J GM* / R* as a function of polytropic 
index, n p i. These two parameters contain the dependence of the surface velocity amplitude 
on the stellar structure. The p-modes have predominantly radial displacement, while g- 
modes have predominantly transverse displacement. 

FIG. 7. — J n 2 as a function of polytropic index, n p i. Rotational splitting of the eigenfre- 
quencies is, to first order in fi*, proportional to J n t (see eq. [16]). For p-modes, J n 2 
approaches as the order n increases. For (/-modes, J n 2 approaches 1/6 as n increases. 

FIG. 8. — Fractional luminosity variation SF"^/ ' F for several low-order (/-modes as a func- 
tion of orbital period T or &. The secondary star is taken to be a n p i = 3 polytrope with 
M* = 5M©, and the primary star has M p = 2Mq. For all the modes considered, 8r a is 
fixed artificially at 0.1, to avoid confusion due to close resonances, and d a = 0.001. The 
eccentricity of the orbit is taken to be 0.8 and 0.4 (Fig. 8a and Fig 8b, respectively). 
The lower limits for T or & of approximately 18 days (for e = 0.8) and 4 days (for e = 0.4) 
correspond to the onset of the overflow of the secondary's Roche lobe. For larger T or &, 
the luminosity variation due to the higher overtone g modes dominates because of their 
smaller b a values. Morever, the variations are significantly larger for the more eccentric 
orbit. 
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